Lamb shift calculated by simple noncovariant method 
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The Lamb Shift (LS) of Hydrogenlike atom is evaluated by a simple method of quantum elec- 
trodynamics in noncovariant form, based on the relativistic stationary Schrodinger equation. An 
induced term proportional to ~p* 4 in the effective Hamiltonian is emphasized. Perturbative calcula- 
tion of second order leads to the LS of 1S 1 i 2 state and that of 2 l S' 1 / 2 — 2P 1 / 2 states in H atom with 
the high accuracy within 0.1% 
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The experimental discovery of Lamb Shift (LS) in 1947 and its theoretical explanations that followed are of great 
importance for the establishment of Quantum Electrodynamics (QED). (Ref. 0- ||). The experimental value of the 
energy difference between 2S\/ 2 and 1P\i-x states for Hydrogen atom reads (in unit of microwave frequency) 

E{2S 1/2 ) - E(2P 1/2 ) = 1057 M5M Hz (1) 



^ ■ while the absolute LS for lSjya state is Q 

AE H (lS 1/2 ) = 8172MMHz (2) 
AE D (lS 1/2 ) = 8184.00M Hz (3) 

G\ 

The LS only accounts for the order of 10~ 6 or 0(a 3 ) of that of the binding energy for electron, i.e., that of Rydberg 
^ ' energy 

^ 1 

R y =R H = -a 2 fx = 3.28805128 x 10 9 MHz (4) 

Ch ' 2 

^ : 

, where \i is the reduced mass of electron and a = ^- = 1/137.0359895. (U = c = 1) 
The theoretical investigation over 50 years reveals that: 

(a) The main contribution of LS comes from the difference of radiative correction (i.e., the perturbative energy 
stemming from emitting a virtual photon and then absorbing it) in different states of electron. 

(b) The difference between the wave functions of 5* states and P states is important. The electron in S states has 
more probability to move into the vicinity of nucleus. In other words, it has more high momentum components in the 
momentum representation of wave function for S states. 

In some literature, in the integration of momentum k of virtual photon, the range of k was often divided into 
two regions. For low k from k = to, say, k = am e = am, the binding effect of electron is taken into account 
in noncovariant theory, whereas for high k up from k = am the covariant theory of QED is applied. This kind of 
treatment seems to us is difficult to avoid the double counting in virtual electron states conceptually. Moreover, the 
so-called long wave (i.e., El) approximation was used in the noncovariant theory at low k region as e lkr ~ 1. But 
exp(ikr) ~ exp(iama) ~ exp(i) ~ 1 (a = 1/a/i being the Bohr radius) is also doubtful to be a good approximation. 

We wish to restudy the problem in noncovariant scheme. First of all, the mystery of LS is not only related to the 
small scale of energy shift shown at Eqs.(|l])-(|3|) but also to the high accuracy of the noncovariant calculation based 
on the Stationary Schrodinger Equation (SSE): 

H i> = evb (5) 
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and 



According to the theory of special relativity (SR), the energy of free electron reads 



E=V^? = , + ^-^ + - (8) 

The magnitude ratio of the third term to the second one is approximately (j^j) ~ (f - )- ~ — ~ 0(a 2 ) in a 
Hydrogenlike atom. So at first sight, the relativistic modification to H would be expected to account for an energy 
decrease of e to order of 10 5 MHz. However, in fact, the LS shows an energy increase of S states to order of 10 3 MHz 
only. This is a mystery we should consider first before the LS could be understood in the noncovariant formalism of 
QED. 

II. RELATIVISTIC SSE 

In Refs. 0]- ||, it was argued that the SSE, Eq.(||), is essentially relativistic as long as the eigenvalue e is related 
to the binding energy B as follows: 

5 = A/c 2 [l-(l + -^n (9) 

where M = m + mjy is the total mass of Hydrogenlike atom with rnjy being the mass of nucleus. Obviously, when 
e <C Mc 2 , B ~ —e as expected. 

Eqs.(||)-(0) together with Eq.(||) was derived from the time-dependent Schrodinger equation of two-body (say 
electron-proton) system in combination with a basic symmetry 

Q(-rt,-rp,-t) = x(rt,rp,t) (10) 

The explanation is as follows. The electron and proton are all not pure. They not only have a particle state 
9 {ft, rj, t), but also have a hiding antiparticle state x{fe,fp,t), 9 and x are coupled together according to the sym- 
metry (|lCj). For one body system, this symmetry leads to Klein-Gordon equation (without spin) or Dirac equation 
(with spin) (Refs. ja]) . The symmetry that "the space-time inversion is equivalent to particle-antiparticle trans- 
formation" shown as ( jl0|) is the essence of SR. The various strange effects of SR are nothing but the reflections of 
antimatter which is in a subordinate status (|%| < \9\) and is just displaying its presence tenaciously |J. 

Therefore, at the level of quantum mechanics, the Hamiltonian Hq in SSE Eq.(||) is enough in the form of Eq.(||). 
No p 4 term like that in Eq.(^) is needed for relativistic correction. The latter is brought into consideration implicitly 
at the final stage e — > — B as shown by Eq.(||). 

III. SEMIEMPIRICAL CALCULATION OF LS 

Based on SSE, we consider that in an effective Hamiltonian a term bip will be induced at QED level by radiative 
correction. It will lead to an energy shift of Hydrogcnlike atomic level as: 



AE%% = (Znl \b 2 p*\ Znl) = [-^- - 3]^_ ( a = J_) (n) 

(21 + 1) nrar a/j, 



The coefficient 62 can be fixed from experiment by the LS between 25 l 1 / 2 — %Pi/2 states in Hydrogen: 



7 q 

-j- = -(1057.845 - 0.087)M Hz = 1586.637MiJz (12) 



'H 
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where 0.087A/ Hz is the small correction stemming from the finite nucleus radius. 

To calculate the absolute LS of 15 state, Eq.(||), we should add three extra contributions: 
(a) The vacuum polarization in QED induces a decrease in charge of electron: 110] , [pTj 



Aa = -—^ 13 



which leads to an increase of energy 



For Hydrogen IS state, it reads 



AE^ P 1S = 271.U0MHz (15) 

(b) The relativistic correction from Eq.(^) reads 

AE --2M-~Z a M^ (16) 

which yields a decrease in energy: 

ae h% = -23.8UMHz (17) 

(c) The correction from finite nucleus radius reads 



AET = 5^ ( T } Ry5m (18) 
which contributes a small increase of energy 

AE%J S = 0.697MHz (19) 

Altogether, we obtain theoretically 

AE H h is rV = 8181.208MHz (20) 
The deviation between ( pp| ) and the experimental value, Eq.(|J) is only 0.1%. 

IV. CALCULATION OF B 2 FROM THE FIRST PRINCIPLE 

We are now in a position to derive the value of coefficient 6 2 from the first principle of QED in noncovariant form. 
Consider an electron with charge — e is moving with (three dimensional) momentum ~]f in the center of mass system 
of Hydrogenlike atom and is coupled to the electromagnetic field via two kinds of interactions ( 0] , Q ) 

#W _ JL jjt • f (21) 

)1C 



ij(2) = . ^ X 1 (22) 

2fic 

According to the perturbation theory in quantum mechanics, the energy shift due to (^) from original e v = 4^- 
will be 
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where the stationary state \~f?) = -j=e l P' r is normalized in a volume V. The intermediate (virtual) state \i) is 
composed of a plane wave eigenstate of SSE (Eq.(^|)) if and a virtual photon with continuous momentum k , see 



Fig. 1. So the ei in the denominator of ( p3| ) reads Si — -|— - + u>k, {<-Ok 
electromagnetic field reads as usual: 



k). The quantized potential A of 



J (2^)3/ 2 V2^T A tT 2 fc ' A V k k 



(24) 



After the integration of matrix element of with respect to space, we substitute one S function S (^jf — if — k 



by (2 V \i ancl then perform the integration with respect to if, yielding 



A £(D = _«p! f 1 dr ,(l-r, 2 )I 
J-i 



p 



(25) 



/ = 



dk 
k + i 



(26) 



where rj = ^ ^ and £ = 2(fi — prf). 



V. THE RENORMALIZATION IS A PROCEDURE TO RECONFIRM THE MASS. 



As in the calculation of QED in covariant form, we also encounter the divergent integral /, Eq.(|26]), here in the 
integration of three dimensional momentum k of virtual photon. To treat the divergence, we will follow the spirit of 
a simple but effective method used in covariant quantum field theory, which evolved from the so-called differential 
renormalization in the literature [l3| - pUjl, then was proposed by Ji-feng Yang [ fL7| and applied extensively in Refs. 
|lg|j , (see also the discussion in |20||). Here the trick is as follows. 

Take the derivative of integral with respect to the parameter £ having a mass dimension: 



81 



dk 



o {k + ir 



which is convergent now. Then we reintegrate Eq.(27) with respect to £ for returning back to I 

J = -ln£ + Ci 



where an arbitrary constant C\ appears. Substituting (|28| ) into (25), one obtains 



-iG(-y + 

7T 3 jl jJL 



av„2,p, ai p ^ )Hl + P_ ) + { l { V_ f _P_ + 1L 
3p [i 3 (1 /i 3p 

9 (i 3 3 3 3 /i 



)ln( 



(27) 



(28) 



(29) 



a 

TTfl 



In 2 



In// — —C\ 



(30) 



7T/i J 



2 
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(31) 



Note that, however, the term b^p 2 will be combined into the kinetic energy term in original Hq. They are 
indistinguishable. The appearance of arbitrary constant C\ precisely reflects the fact that we can not calculate the 
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reduced mass of an electron via the perturbation calculation of AEp 1 ^ . So the value of C\ is chosen such that = 0, 
implying that the value of reduced mass /i (yet not /J, b s , see below) is reconfirmed as an observed mass which can 
only be fixed by the experiment (not by theory). 

Next turn to H^ 2 \ which induces the spin flip process between |??,±i) and |lf,±2) states, yielding 



\(i\HW\T?,s z )\ 



dijj 



(32) 



J = 



oo ,2 



k 2 dk 



For this divergent integral, derivative of third order is needed to render it convergent: 

<9 3 J _ 2 

So after reintegrating with respect to £, we have: 

J = -i 2 m£ + C 2 £ 2 + C 3 £ + C 4 



(33) 



(34) 



(35) 



= —{^[(1 + — ) 3 ln(l + — ) — (1 — — ) 3 ln( 
' 7T 3p n [1 ii 



l- P - 



)\--i P -f-- 
n 9 V 3 



+4(ln2 + ln M )-4C 2 -^-% + f^ln2 + 2 + ^ln Ai -^C 2 ") (^) 2 } 



(36) 



(2) _ a/i 

7T 



9C C 
4(ln2 + ln M )-4C 2 - — 



(37) 



&< 2 > = " 



4 4 4 

-ln2 + 2+-ln M --C 2 
TTfi \3 3 3 



(38) 



,(2) 
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(39) 



We shall fix three arbitrary constants C 2 , C3 and C4 carefully. First look at the b^p 2 term which should be 

2 

combined with the term |^ with /i already fixed. Any modification on [i must be finite and fixed. So the only possible 
choice of C 2 is to cancel the ambiguous term | ln,u in b^\ leaving 



1 2/x 



2a /4 



-Zn2 + 2 



(40) 



(2) 

The constants C3 and C4 must be chosen such that 6J, ; 
point. There is always no rest energy term in SSE. 



0, which means that we reconfirm the SSE as our starting 



(2 

Hence, the nonzero contribution of b\ ! p 2 docs bring a finite and fixed modification on fi so that 



i^obs 



l + (3 



(41) 
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where /x b s is the reduced mass eventually observed in experiment. However, as discussed in Eq.(||), there is no term like 

4 

— „ p a in SSE. The rclativistic correction is left to the modification of e to B at the final stage. So an unobservable 

term is also modified, the difference — g (^t -j^j P 4 should be treated also as an unobservable background of 

modification in p 4 term, which should be subtracted off from the observable one. Hence while b\ = 6j + = b] 2 , 
but instead of 62 = &2 + ^ , 2 2 ' > = ^3 ( — i)i we should have a renormalized b R as 

& * = & 2 + s-s (3/3 + 3/3 2 + /3 3 ) = — V (1-942816878) (42) 



VI. COMPARISON WITH THE EXPERIMENTAL VALUES 

Let us use Eq. to evaluate the LS of H atom: 

Kh R 

*Eb& = {b2P 4 ) ls = -r- = 7901.629MH* (43) 



Adding the contributions from Eqs.(|l5|), (fi7|) and fll9|), we get 



AEl h ™ ry = 8149.653MHz (44) 



which is smaller than the experimental value Eq.(||) up to 0.28%. 
For the LS between 2S , 1 / 2 and 2P\/2 states, 



2 b R 

^E R % 2P = --£- = 1053.551MHz (45) 



After adding a small contribution due to the finite nucleus radius, we have 



AE H h 2s r J 2 -2P 1/2 = 1053.638MHz (46) 



which is also smaller than the experimental value Eq.([j]) up to 0.40%. 

For further improvement, we keep all p n (n > 4) terms. So we manage to evaluate the renormalized radiation 
correction as follows 

AE Rad {~? ) = A£« + AE& _ 6< 2 V _ (Jp2 +// 2 _ ^ obs _ JL) + (v ^T^ - H - (47) 



with /^ h s — = d[i = —/3^ bs and \/p 2 + n 2 obs — \/p 2 + [J. 2 = ^yMf as we wish to keep the explicit dependence on 

a throughout Eq. ( |47| ) being of order O(a). Then we calculate the expectation value of AE Rad (J?) for a fixed state 
numerically, yielding 

AE§^ S = 7920.533MHz (48) 

A^2S 1/2 -2P 1/2 = 1057.550MHz (49) 

A£g^ s = 7922.688MHz (50) 

A££ a 25 1/2 - 2 p 1/2 = 1057.838MHz (51) 
After adding the other three corrections Eqs.(p^|)-(p^|), we are pleased to see the results 

AE l h is ry = 8168.557MHz (52) 
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AE l h 2sJ 2 -2P 1/2 = W57M7MHz (53) 

AE l h is rV = 8186.181MHz (54) 

Ae d!2sI V 2 -2p 1/2 = 1058.363MHz (55) 
coinciding with the experimental data to a high accuracy (< 0.1%). 

VII. SUMMARY AND DISCUSSION 

(a) We propose a simple but effective method for calculating the LS in a noncovariant form. Two crucial observations 
are as follows: 

(i) The SSE is essentially relativistic as long as its eigenvalue e is related to the binding energy B by Eq.(^). 

(ii) The main contribution to LS is coming from the different mean values of operator "J? 4 in S and P states as shown 
in the semiempirical calculation. The analysis convinced ourselves that the use of three dimensional momentum J? is 
much suitable than that of four dimensional one. Actually, a simple one-loop calculation in covariant form of QED 
gave us the value of LS for H atom with accuracy only 5% for 2£ , 1 / 2 — states and even worse (~ 22%) for lS'1/2 
state 0. 

(b) The subtlety of SSE can be seen further by the following comparison. At free moving condition, the relativistic 

2 

effect is contained in the definition e = ^ with rest (reduced) mass being a constant containing all the radiative 
correction. On the other hand, when the electron is bound in a Hydrogenlike atom, its Hamilton H$ contains no rest 
mass and no explicit (negative) 7? term either. The relativistic effect is contained in the definition e = — ^M- (with 
E and M being the total energy and mass of the system) and Eq.(j^) implicitly. However, the radiative correction 
does induce a small (positive) term of ~jf A at the level of QED. 

(c) The definition of s and Eq.(^|) also indicate that there is no any negative energy eigenstate in SSE. The hiding 
antiparticle state (which is the essence of special relativity) is already taken into account in deriving SSE with Eq. ([)]). 
In other words, no further explicit virtual positron state should be considered in our calculation. Actually, we had 
struggled for years before eventually realizing that only the simple formulas ( p3| ) and ( |32|) with Fig.l are needed in 
the perturbative calculation of second order. 

(d) The effective Hamiltonian of Hydrogenlike atom for evaluating the LS can be summarized as 

«■><= <» 

while the third term (~ b^p 4 ) is the QED modification to the first term, the fourth term could be seen as that to the 
second term (see Eqs.([l3|), 0)). 

Eq.(p6|) is applicable to j = 1 states. The P3/2 state is pushed up by extra spin-orbit coupling to form the fine 
structure. For Hydrogen, E(2P 3/2 ) - E(2P 1/2 ) = 1.09691 x 10 4 MHz about 10 times of LS. Furthermore, the hyperfine 
structure (hfs) in Hydrogen, stemming from the interaction between the magnetic moments of electron and proton, 
is smaller than the LS. The energy splitting of lS'1/2 states is well known as AEis 1/2 (hfs) — 1420.406AfiJz. As we 
don't take the magnetic moment of nuclei into account, the hfs is not considered in this paper. 

In deriving the H e f / , the Pauli interaction term H ^ between electron spin and the external magnetic field has to 
be added to the RSSE so that the spin-orbital coupling (hyperfine structure) and LS can be calculated quantitatively. 
This is the price we must pay for the use of RSSE. On the other hand, though the Dirac equation in external field can 
predict the electron spin with g = 2, it fails to take the difference of reduced mass of electron into account because it 
is a one-body equation. Furthermore, it predicts a too low ground state (IS) and a too large splitting of 2P 3 / 2 and 
2-Pi/2 states. The weakness of Dirac equation seems to us is due to its overestimation of the antiparticle ingredient 
in the electron as discussed in Ref. Actually, only after long time study on Dirac equation with its difficulty 
in calculating LS, especially the absolute LS of IS states, could we believe in the advantage of using RSSE as the 
starting point. 

(e) The reasonable result in this paper again shows the correctness and effectiveness of the renormalization method 
used in Refs. |l7]]- pQ] . For treating the divergence , which warns us of the lack of our knowledge about the parameters 
(mass, charge, etc.), the renormalization is a procedure to reconfirm the parameters step by step rigorously. Here 
it is interesting to see an example of finite and fixed renormalization. Note that, however, we must consider the 
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contribution of first to define the parameter jj, before to consider that of for bringing [i to fi bs- This is the 
only reasonable logic. The inverse logic, i.e., to consider first and then next, would lead to inconsistent and 
wrong result. 
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FIG. 1. The only one loop diagram calculated for radiation correction in the noncovariant formalism of QED. At the two 
vertices, either or is used (no interference between them would occur due to no polarization in the plane wave p 

state) . 
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